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Abstract 

We calculate the branching ratios and CP-violating asymmetries for vr'^r/^'^ decays in the 

perturbative QCD (pQCD) factorization approach here. We not only calculate the usual factor- 
izable contributions, but also evaluate the non-factorizable and annihilation type contributions. 
The pQCD predictions for the CP-averaged branching ratios are BR{Bg — > tt^t]) « 0.86 x 10"'^ 
and BR{B^ — > -n^rj') ~ 1.86 x 10~^. The pQCD predictions for the CP-violating asymmetries 

are A'^^'p{-K^y]) 4.5%, ^^^^(vrO??') 9.1%, A'gfi^^r]) 0.2%, and A'gf{-K^'n') ~ 27.0% 

but with large errors. The above pQCD predictions can be tested in the near future LHC-b 
experiments at CERN and the BTeV experiments at Fermilab. 
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I. INTRODUCTION 



The experimental measurements and theoretical studies of the two body charmless 
hadronic B meson decays play an important role in the precision test of the standard 
mode (SM) and in searching for the new physics beyond the SM P|. For these charmless 
B meson decays, the dominant theoretical error comes from the large uncertainty in 
evaluating the hadronic matrix elements {MiM2\0i\B) where Mi and M2 are light final 
state mesons. The QCD factorization (QCDF) approach and the perturbative QCD 
(pQCD) factorization approach 0,1^ are the popular methods being used to calculate the 
hadronic matrix elements. 

When the LHC experiment is approaching, the studies about the decays of me- 
son draw much more attentions then ever before. At present, some two-body charmless 
hadronic Eg meson decays have been calculated, for example, in both the QCDF ap- 
proach and/or in the pQCD approach In this paper, we would like to calculate the 
branching ratios and CP asymmetries for Bg — > tt^t^^'^ decays by employing the low energy 
effective Hamiltonian and the pQCD factorization approach. Besides the usual factor- 
izable contributions, we here are able to evaluate the non-factorizable and the annihilation 
contributions to these decays. 

Theoretically, the two Eg — > tt'^?]'^') decays have been studied in the naive and gener- 
alized factorization approach |^ or in the QCD factorization approach ^| . On the 
experimental side, only the poor upper limits for the branching ratios are available now 
3 

BR{Bl 7r°r/W) < 1.0 x 10"^ (1) 

Of course, this situation will be improved rapidly when LHCb starts to run at the year 
of 2007. 

For Eg n^r]^'^ decays, the light final state mesons are moving very fast in the rest 
frame of Eg meson. In this case, the short distance hard process dominates the decay 
amplitude, while the soft final state interaction is not important for such decays, since 
there is not enough time for light mesons to exchange soft gluons. Therefore, it makes the 
pQCD reliable in calculating the Eg — > tt'^t]^'^ decays. With the Sudakov resummation, 
we can include the leading double logarithms for all loop diagrams, in association with 
the soft contribution. 

This paper is organized as follows. In Sec. m we calculate analytically the related 
Feynman diagrams and present the various decay amplitudes for the studied decay modes. 
In Sec. IIIH we show the numerical results for the CP-averaged branching ratios and 
CP asymmetries of Eg n^r]^ ^ decays and compare them with the measured values or 
the theoretical predictions in QCDF approach. The summary and some discussions are 
included in the final section. 
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II. PERTURBATIVE CALCULATIONS 



For Bg n T]^'' decays, the related weak effective Hamiltonian H^jf can be written as 

i 



^eff 



Gp 

71 
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i=3 



(2) 



The exphcit expressions of the operators Oj can be found for example in Refs. |12l . lie 
In the pQCD approach, the decay amplitude is conceptually written as the convolution, 

A{B, ^ M1M2) ~ Jd%d%d% TT[CmBAki)^MAk2)<^M,{h)H{h,k2,ks,t)], (3) 

where fcj's are momenta of light quarks included in each mesons, and Tr denotes the 
trace over Dirac and color indices. C{t) is the Wilson coefficient which results from the 
radiative corrections at short distance. The function H{ki, k2, k^jt) describes the four 
quark operator and the spectator quark connected by a hard gluon whose is in the 
order of AMb^, and includes the 0{^/AMb^) hard dynamics. Therefore, this hard part H 
can be perturbatively calculated. The function is the wave function which describes 
hadronization of the quark and anti-quark to the meson M. While the function H depends 
on the processes considered, the wave function $m is independent of the specific processes. 
Using the wave functions determined from other well measured processes, one can make 
quantitative predictions here. 

Since the b quark is rather heavy we consider the Bg meson at rest for simplicity. It is 
convenient to use light-cone coordinate {p~^,P~,Pt) to describe the meson's momenta, 

= ±p'), and pT = {p\p'). (4) 

Using these coordinates the Bs meson and the two final state meson momenta can be 
written as 

Pi = ^(1,1,0t), P2 = ^(1,0,Ot), P3 = ^(0,1,Ot), (5) 

respectively, here the light meson masses have been neglected. Putting the light (anti-) 
quark momenta in P^, 7r° and t]"-'^ mesons as fci, ^2, and ^3, respectively, we can choose 

fci = (a;iP+,0,kiT), ^2 = (x2P2+,0,k2T), fc3 = (0,a;3P3-,k3T). (6) 
Then, the integration over /c]", A;^, and k'^ in eq.© will lead to 



A{Bs ^ n'^r]'^'^) ~ I dxidx2dx^hidhih2dh2h2,db 



'3 



■Tr 6i)$.o(x2, 62) V) (^3, h)H{x,, b„ t)Stix,) e~^^'>] (7) 

where bi is the conjugate space coordinate of kiT, and t is the largest energy scale in 
function H{xi,bi,t). The large double logarithms (In^Xj) on the longitudinal direction 
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are summed by the threshold resummation and they lead to St{xi) which smears 
the end-point singularities on Xi. The last term, e"'^'-*-*, is the Sudakov form factor which 
suppresses the soft dynamics effectively [l^. In numerical calculations, we use = 

47r/[/5i \n{t'^ /AqIjj^ )] which is the leading order expression with AgJ,^ = 193MeV, derived 

from Aq^^ = 250MeV. Here /3i = (33 — 2nf)/12, with the appropriate number of active 
quarks n/. 

Similar to B ^ pr]^'^ and B tct]^'^ decays, there are 8 type diagrams contributing to 
the Bs — >■ TT^T]*^'^ decays, as illustrated in Figure 1. We first calculate the usual factorizable 
diagrams (a) and (b). Operators Oi, O2, O3, O4, O9, and Oio are (V — A)(y — A) currents, 
the sum of their amplitudes is given as 

/• 1 /"OO 

Ferj = 4:V27TGFCFfnm%^ / dXidX3 / 61^6163^63 </)B^ (Xi, 61) 

Jo Jo 

X { [(1 + X3)0^ (X3, 63) + (1 - 2x3)r^(0^(x3, 63) + 0J^(X3, 63))] 

■as{tl) heixi,X3, 61, 63) exp[-Sabitl)] 
+2r'^(j)^{x3, h3)as{tl)he{x3, xi, 63, 61) ey.^[-Sab{tl)] } , (8) 

where = itiq"" /ttib/, Cp = 4/3 is a color factor. The function he, the scales and the 
Sudakov factors Sab are displayed in AppendixEl In the above equation, we do not include 
the Wilson coefficients of the corresponding operators, which are process dependent. They 
will be shown later for different decay channels. 

The form factors of Bg to t]*-'-* decay, F^f^^ (0), can thus be extracted from Eq. (jH)), 
that is 

F^r'"\q' = 0) = — • (9) 

The operators O5, Oe, O7, and Os have a structure of {V — A){y + A). In some decay 
channels, some of these operators contribute to the decay amplitude in a factorizable way. 
Since only the axial- vector part of {V + A) current contribute to the pseudo-sealer meson 
production, {'k\V - A\B){'q\V + A|0) = -{■k\V - A\B){'q\V - A|0), that is 

F^n = -Pen ■ (10) 

For the non-factorizable diagrams 1(c) and 1(d), all three meson wave functions are 
involved. The integration of 63 can be performed using S function 5(63 — 62), leaving only 
integration of 61 and 62. For the (V" — A)(y — A) operators, the result is 

16 /'CO 

= —='nGFCFrrL^Bs / dxidx2dx3 / 61^6162^62 0b^(xi, 61) 
v3 Jo Jo 

X0^(X2, 62)^3 [(p^{x3, 62) - 2r^0j'(x3, 62)] 

■asitf)hf{xi,X2,X3,bi,b2)exp[-Scd{tf)] ■ (11) 

M^^ is for the {S — P){S + P) type operators, which are from Fierz transformation for 
{V - A){V + A) operators: 

< = -M,,. (12) 
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For the non-factorizable annihilation diagrams 1(e) and 1(f), again all three wave 
functions are involved. Here we have two kinds of contributions. M^r, and describe 
the contributions from the {V — A)(y — A) and {S — P){S+P) type operators, respectively, 

16 /'OO 

Mari = -^-nGpCFrn^Bs / dxidx2dx3 / bidbib2db2 (f)Bsixi^bi) 

v3 Jo Jo 

X { {X3(P^{X3, b2)(pt{x2, &2) + V^''^ [x2 (^f (xa, ^2) - ^1{X2, 62)) " 

(0j'(x3, 62) - 0j^(a;3, ^2)) + 3:3 (0f (a;2, ^2) + 0^(a;2, ^2)) ■ 
{(f)^{x3, 62) + 0^(3^3, ^2))] } as(t])h]-{xi, X2, X3, 61, 62) exp[-S'e/(t})] 
- {x20;f (x3, &2)0;^(a;2, &2) + V^;-'^ [((a;2 + 0:3 + 2)0j'(x2, 62) + (^2 - X3) 
x0j(a;2, 62)) 0^(a;3, ^2) + {{x2 - X3)0^(x3, 62) + {^2 + X3 - 2) 
0j(a;3,&2)) 0^(a;2, 62)] } as(i/)/i/(a;i, 0:2, X3, 61,62) exp[-S'e/(tJ)]} , (13) 
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16 /"CXD 

= —j=TiGFCF'rrt^Bs / dxidx2dxz I hidhih2dh2(l)Bs{^ii^i) 
v3 Jo io 

X { {X20;^(X3, &2)0;f (^2, &2) + V^''^ [xg (0^(X2, 62) - 0^X2, 62)) " 

«(a;3, &2) - <(X3, &2)) + a;2 (^f (0:2, 62) + 0na;2, &2)) ■ 
(0^(x3, 62) + 0^(3:3, 62))] } as{t))h){xi, X2, X3, bi, 62) exp[-S'e/(t})] 
- {3:20^(0:3, &2)0;^(X2, &2) + r.r^''^ [((X2 + 0:3 + 2)0J'(X2, 62) + (X3 - X2) 

X0J'(X2, &2)) 0f (a;3, &2) + ((3^3 " X2)(p^{Xs, 62) + {^2 + X3 - 2) 

0^(3:3, ^2)) (Pl{x2, &2)] } as(t5)/ij(a;i, 0:2, X3, 61, 62) exp[-S'e/(tJ)]} . (14) 

where r;;''^' = mo"^''7mB,. 

The factorizable annihilation diagrams 1(g) and 1(h) involve only 7r° and r]^'^ wave 
functions. There are also two kinds of decay amplitudes for these two diagrams. is 
for (y — A)(y — A) type operators, is for (V — A){V + A) type operators: 

/•l /"OO 

= Far, = 4:V2'KGFCFfBs'm%l dX2dX3 / b2db2hdb3 {[X3(f)^{x3,b3)(f)^{x2,b2) 

Jo Jo 

+2v!;^\{x3 + 1)0J'(X3, &3) + (3^3 - l)0j(a:3, &3))0f (a:2, &2)] 
■as{tl)ha{x2, X3, 62, &3) exp[-S'g/i(tg)] 

- [2^20^(a;3,&3)0;^(a;2,&2) 

+2r^r;''^0j'(x3,63)((x2 + 1)0^(^2,62) + (^2 - 1)0^2^2,62))] 
•a^(tg)/ia(x3,a;2,63,62)exp[-S'gh(tg)]} . (15) 

If we exchange the vr and r]^'^ in Fig. 1, the corresponding expressions of amplitudes 
for new diagrams will be similar with those as given in Eqs.(jSllT3j). The expressions of 
amplitudes for new diagrams can be obtained by the replacements. 



0^-0j', 0^-0j, r.'-^'" 



r;''. (16) 



For example, we find that 

Fa7T = -Far^C), = -F^riC) ■ (17) 

Now we are able to calculate perturbatively the form factors Fq "'^^ (0) and the decay 
amplitudes for the Feynman diagrams after the integration over Xi and 6,. Since we here 
calculated the form factors and amplitudes at the leading order ( one order of dsit)), the 
radiative corrections at the next order would emerge in terms of as{t) ln(m/t), where m's 
denote some scales, like niB^, l/h, . . ., in the hard part H(t). We select the largest energy 
scale among m's appearing in each diagram as the hard scale t's for the purpose of at 
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least killing the large logarithmic corrections partially, 

tl = ■ max(A/i^mij^, 1/61, 1/63) , 
tl = at- max(yxrmB,, I/&1, 1/63) , 
tl = at-max^y/x^mB^l/hA/h) , 



= at■max{^/x^mB,,l/b2,l/b. 



'3j 



at ■ max{y/x^mB,, \/x^2,mB,, l/hi, I/62) , 
t) = at ■ msix{y/x2X3mB,, l/h, I/62) , 

tj = at ■ max{y/xi + X2 + X3 - XiX^ - X2X3mB,, y/x2X3mB,, 1/bi, 1/62) , (18) 

where the constant at = 1.0 ± 0.2 is introduced in order to estimate the scale dependence 
of the theoretical predictions for the observables. 

In Ref.jl^ [l3| . a brief discussion about the 77 — 77' mixing and the gluonic component 
of the 77' meson have been given. Here we don't show it again. 

Combining the contributions from different diagrams, the total decay amplitude for 
TT^f] can be written as 

V6Min%) = Fe, (c, + ^^2^ - 6 (-^Cr - + ^C, + ^C,,^ | ^2(0,) 

|e„C2 - 6 (-^Q + lc^o^ } F2{9,) + (^M,, + 

aC2 - e^C^io} F,ie,) - 6 (m^; + M;^^ ^C8Fi(^,) . (19) 

The decay amplitudes for —>■ n^r]' can be obtained easily from Eqs. lfT^ by the 
following replacements 

sin ^ 



V2 ' 

F2i9p) — > F^{9p) = cos 9p- V2 sin 9p. (20) 
Note that the possible gluonic component of rj' meson has been neglected here. 



III. NUMERICAL RESULTS AND DISCUSSIONS 
A. Input parameters and wave functions 

We use the following input parameters in the numerical calculations 

Agr^) = 250MeV, ^ = 130MeV, /b, = 230MeV, 
m'^o"'' = 1.4GeV, ml'' = 1.95GeV, fx = 160MeV, 
Mb, = 5.37GeV, Mw = 80.41GeV. (21) 

For the CKM matrix elements, here we adopt the Wolfenstein parametrization for the 
CKM matrix, and take A = 0.22, A = 0.853, p = 0.20 and r] = 0.33 HH]. 
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For the Bs meson wave function, we adopt the model 



Ml 



X 



2-1 



MM' 



(22) 



where is a free parameter and we take = 0.50 ±0.05 GeV in numerical calculations, 
and Nb^ = 63.7 is the normalization factor for uJbs = 0.50. 

For the light meson wave function, we neglect the b dependant part, which is not 
important in numerical analysis. We use the wave functions of vr meson ( 0^(a:), (p^ix) 
and (f)'^{x) ) as given in Ref.jlq]. For rj meson's wave function, and represent 

the axial vector, pseudoscalar and tensor components of the wave function respectively, for 
which we utilize the result from the light-cone sum rule including twist-3 contribution. 
For the explicit expressions of the wave functions and the values of related quantities, one 
can see Eqs.(50) and (51) of Ref. . 

We assume that the wave function of uu is same as the wave function of dd. For the 
wave function of the ss components, we also use the same form as dd but with an d fy 
instead of ttiq'^ and fx, respectively. For and fy, we use the values as given in Ref.nJ 
where isospin symmetry is assumed for fx and SU{3) breaking effect is included for fy-. 



fx fjT, fy 



(23) 



These values are translated to the values in the two mixing angle method, which is often 
used in vacuum saturation approach as: 



/s = 169MeV, A 
^8 = -25.9°(-18.9'= 



151MeV, 

= _7.r(-o.r). 



(24) 



where the pseudoscalar mixing angle Op is taken as —17° (—10°) [19]. The parameters m, 
= Vddiuu),Vss) are defined as: 



m, 







rrin 



{ruu + ma) 



2Ml - ml 
(2m,) 



(25) 



We include full expression of twist— 3 wave functions for light mesons. The twist— 3 
wave functions are also adopted from QCD sum rule calculations jl^. We will see later 
that this set of parameters will give good results for Bs tt^t]^'^ decays. 



B. Branching ratios 

For Bs — * TT^T]^'^ decays, the decay amplitudes in Eqs. (|T!H) can be rewritten as 
M = V:,VusT - V:,VtsP = V:,VusT [l + ze^(^+^)] , 



where 



(26) 



VtlVts 




p 


v:bVus 




T 



(27) 
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is the ratio of penguin to tree contributions, 7 = arg — is the weak phase (one of 

us 



the three CKM angles), and 6 is the relative strong phase between penguin (P) and tree 
(T) diagrams. In the pQCD approach, it is easy to calculate the ratio z and the strong 
phase 6 for the decay in study. For Bg —>■ tc^t] and tt^i]' decays, we find numerically that 

z(7r%) = 38.3, 5(A) = -94°, (28) 
zi-K^T]') = 5.5, 5(7rV) = -20°. (29) 

The main error of the ratio z and the strong phase S is induced by the uncertainty 
of = 0.50 ± 0.05 GeV. Since the errors induced by the uncertainties of most input 
parameters are largely canceled in the ratio, we will use the central values of z and 6 in 
the following numerical calculations, unless explicitly stated otherwise. 

Using the wave functions and the input parameters as specified in previous sections, it 
is straightforward to calculate the branching ratios for the four considered decays. The 
theoretical predictions in the pQCD approach for the CP-averaged branching ratios of 
the decays under consideration are the following 

Bri 5° ^ vr%) = [0.86l°iI(a;,J+°i?(m.)+J:°°(a,)] x 10"^ (30) 
Bri 5° ^ vr^O = [l.86l°f (a;,Jl°f x 10-^ (31) 

for Op = -17°, and 

Br{ Bl 7r%) = [l.l8l°:i(o;,Jl°:i(m.)lJ:?^(a,)] x 10"^ (32) 
Br{ 5° ^ vr^O = [l.54;°;^^(o;,J+°:i(m.)+^i?(a,)] x 10"^. (33) 

for 9p = —10°. The main errors are induced by the uncertainties of at = 1.0 ± 0.2, 
Ub, = 0.50 ± 0.05 GeV and m, = 120 ± 20 MeV, respectively. 

It is easy to see that (a) the errors of the branching ratios induced by varying at in the 
range of at = [0.8, 1.2] can be significant for the penguin-dominated B^ n'^T]^'^ decays; 
and (b) the variations with respect to the central values are large for the case of at = 0.8, 
but very small for the case of at = 1.2). This feature agrees with general expectations: 
when the scale t become smaller, the reliability of the perturbative calculation of the form 
factors in pQCD approach will become weak! 

The pQCD predictions of the branching ratios as given in Eqs. ()30llS^ agree well with 
the theoretical predictions in the QCDF approach, for example, as given in Ref. 

Br{ 5° ^ A) = (0.75^°:i) X 10-^ 

fir(fi°^A') = (1.11°:^^) X 10-^ (34) 



where the individual errors as given in Ref. [5[ have been added in quadrature. 
C. CP- violating asymmetries 

Now we turn to study the CP-violating asymmetries for 5^ -k^t]^'^ decays. For these 
neutral decay modes, the effects of 5° — B^ mixing should be considered. 
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Acp = 



For Bg meson decays, we know that AT/Amg -C 1 and AT/T <^ 1. The CP-violating 
asymmetry of 5^(5^) — > tc^t]^'^ decay is time dependent and can be defined as 

r (50(At) ^ fcp) - r (i?o(At) ^ fcp) 

r (BO{At) ^ /cp) + r {BO{At) ^ /cp) 

= A'^jp cos(Am,At) + A^)? sin(Am,At), (35) 

where Am^ is the mass difference between the two mass eigenstates, At = tcp — ttag 
is the time difference between the tagged 5° {Bg ) and the accompanying Bs (5°) with 
opposite b flavor decaying to the flnal CP-eigenstate fcp at the time tcp. The direct and 
mixing induced CP-violating asymmetries and can be written as 

Adir _ \^cp? - 1 .mix _ '^Im{Xcp) , . 

-^^^ " l + |AcPr " l + |AcPr 

where the CP- violating parameter Acp is 

Here the ratio z and the strong phase 5 have been deflned previously. In pQCD approach, 
since both z and 5 are calculable, it is easy to flnd the numerical values of A*p and 
for the considered decay processes. 

In Figs. El we show the 7— dependence of the direct CP- violating asymmetry A^*p for 
TT^T] (solid curve) and B^ n^rj' (dotted curve) decays for 9p = —17° . 

The pQCD predictions for the direct CP-violating asymmetries of 5° tx^t]^'^ decays 

are 

A'^p{Bl TT^) = ± ^.Q{ml)t\l{m,ni{a,)\ x (38) 

^^-(5°^ AO = [-9-l-2:^(7)^o:6K)±0.3(m^)±1.9(m.)lli(a,)] xlO-l (39) 

As a comparison, we present the QCDF predictions for ^^*p(i?° tt^t]') directly 
quoted from Ref. 

A'MB's - vr^O = (27.81?:? To tl^'l) x 10-^ (40) 

where the "default values" of the input parameters have been used in Ref. , and the error 
sources are the same as the flrst four input parameters in Eqs. (jHHj) and (jH^ . Currently, 
no relevant experimental measurements for the CP-violating asymmetries of B^ tt'^t]^'^ 
decays are available. For the direct CP-violating asymmetries of B'^ tt'^t]^'^ decays, the 
theoretical predictions in pQCD and QCDF approach have the opposite sign, but the 
theoretical errors are clearly too large to make a meaningful comparison. One has to wait 
for the improvements in both the experimental measurements and the calculation of high 
order contributions. 

The pQCD predictions for the mixing induced CP-violating asymmetries of B^ tt^t]^'^ 
decays are 

^cp(i?°-A) = [-0-2±0.1(7)^l?K)^ii«)lt5K)^n:6M xlO-^ (41) 
^CP(5°-A') = [27.0l^,|(7)l°:^(c.,Jl°j(m^)±0.2(m,)l^y(a,)] xlO-^ (42) 
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FIG. 2: The direct CP asymmetry Af^p (in percentage) of Bg — > tt^i] (solid curve) and Bg Tr^rj' 
(dotted curve) as a function of CKM angle 7 for the case of 6p = —17°. 
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FIG. 3: The mixing induced CP asymmetry A^^p (in percentage) of Bg — * vr^r/ (solid curve) 
and Bg tt^t]' (dotted curve) as a function of CKM angle 7 for the case of = — 17° . 

where the dominant errors come from the variations of Uf,^ — 0.50 ± 0.05 GeV, ttIq = 
1.4 ± 0.3 GeV, at = 1.0 ± 0.2, = 120 ± 20MeV and 7 = 60° ± 20°. 

If we integrate the time variable t, we will get the total CP asymmetry for — > 7r°?7*^') 
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decays, 

where x = Am^/F = 26.5 for the 5^ — Bg^ mixing We found numerically that the 
magnitude of the total CP asymmetry for n^r]^'^) decays are smaller than 2% in 

the whole considered parameter space. 



D. Effects of possible gluonic component of r]' 

Up to now, we have not considered the possible contributions to the branching ratios 
and CP-viola ting asymmetries of 5° — >■ iri]' decays induced by the possible gluonic com- 
ponent of rj' [l9l UM- When Z^/ 7^ 0, a decay amplitude A4' will be produced by the 



gluonic component of rj'. Such decay amplitude may construct or destruct with the ones 
from the qq {q = u, d, s) components of rj', the branching ratios of the decays in question 
may be increased or decreased accordingly. 

Unfortunately, we currently do not know how to calculate this kind of contributions 
reliably. But we can treat it as an theoretical uncertainty. For \M' /M{qq)\ ~ 0.1 — 0.2, 
for example, the resulted uncertainty for the branching ratios as given in Eq. (j31|) will be 
around twenty to thirty percent. 

Furthermore, the pQCD predictions for the branching ratios oi B ^ p?]^'^ and B — 
7r?7*^') decays also show very good agreement with the data jl^. Es^. We therefore believe 
that the gluonic admixture of r]' should be small, and most possibly not as important as 
expected before. 

As for the CP-violating asymmetries of B^ — > tt^t]' decays, the possible contributions 
of the gluonic components of the t]' meson are largely canceled in the ratio. These results 
may be measured in the forthcoming LHCb experiments. 



IV. SUMMARY 

In this paper, we calculate the branching ratios and CP-violating asymmetries of B^ —>■ 
n^r], Bg —>■ 71^1]' decays in the pQCD factorization approach. 

Besides the usual factorizable diagrams, the non-factorizable and annihilation diagrams 
are also calculated analytically. Although the non-factorizable and annihilation contribu- 
tions are sub-leading for the branching ratios of the considered decays, but they are not 
negligible. Furthermore these diagrams provide the necessary strong phase required by a 
non-zero CP-violating asymmetry for the considered decays. 

From our calculations and phenomenological analysis, we found the following results: 

• The pQCD predictions for the form factors are FqI^^'{0) = —0.276 and F^f^^ (0) = 
0.278, which agree well with those obtained from other methods. 

• For the CP-averaged branching ratios of the considered decay modes, the pQCD 
predictions for Op = 17° are 

Br{B', A) = {OMtlfs) X 10-^ 

5r(5° A') = {l-SQtlit) X 10-^ (44) 
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here the various errors as specified in Eqs. and (pTT|) have been added in quadra- 
ture. The pQCD predictions are also well consistent with the results obtained by 
employing the QCD factorization approach. 

• For the CP- violating asymmetries, the pQCD predictions for Aqp{Bs n'^rj^'^) 
and Acp{Bs n^i]^'^) are generally not very large, while the time-integrated CP 
asymmetries are less than 2% in magnitude. 

• The major theoretical errors of the computed observables are induced by the un- 
certainties of the hard energy scale t/s, the parameters and m^, as well as the 
CKM angle 7 for CP asymmetries. 
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APPENDIX A: RELATED FUNCTIONS 

We show here the function hiS, coming from the Fourier transformations of H^^\ 

he{xi,X3, bi, 63) = Kq {^JxiX^niBbi) [d{bi - b^)Ko {^/x^mBbi) Iq {^/x^niBb^i) 

+9{b3 - bi)Ko {^fximsbz) Iq {^/ximBbi)] St^xs), (Al) 

ha{x2, X3, 62, ^3) = Kq [iy/x^niBbs) [9{b3 - b2)KQ (iy/x^rriBbs) Iq {iy/x^mBb2) 

+9{b2 - b3)Ko {iy/x^mBb2) Iq (iy/x^mBba)] St{x3), (A2) 



hf{Xi,X2,X3,bi,b2) = \9{b2 - bl)lQ{MB^/x^ibl)'^oiMB^/x^b2) 

Ko(MbF(i)6i), for F^) > 0^ 
+ ih - 62) M ^H«(M^^7f^ b,), for Ff,) < 1 ' (^3) 



hf{xi, X2, X3, 61, 62) = <^9{bi - b2)KQ{i^/x^blMB)lQii^/x^b2MB) + {pi ^ 62) 



^-H^ (Va^i + X2 + X3- X1X3 - X2X3 &iMij), (A4) 



hf{xi,X2,X3,biM) = <9{bi -b2)Ko{iy/x^blMB)lQ{i^/x^,b2M^ 



B) 



Ko(MbF(2)&i), for F2)>0^ 
+ ih - 6.) M ^hW(m,/^ bi), for < 1 ' (^5) 
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where Jo is the Bessel function and Kq, Iq are modified Bessel functions Kq{ 
— {7r/2)Yo{x) + i{7i/2)Jo{x), and -F(j)'s are defined by 



1^ 

^(2) 



[Xi - X2)X3 
(Xi - X2)X3 



The threshold resummation form factor St{x-i) is adopted from Ref.j2 

2i+2T(3/2 + c) 

where the parameter c = 0.3. This function is normahzed to unity. 
The Sudakov factors used in the text are defined as 



Sab{t) = s {xiitlb/ V2, 6i j + s {x^itlb/ V2, 63) + s (^(1 - X3)mB/ V2, h 



1 

Ti 



In + In . 



-ln(6iA) -ln(63A) 
s {xiItlb/ V2, 61 j + s (^X2mB/ V2, b 

+s {x^mB/ v^, &2 j + s (^(1 - X3)mB/ V2, 62 j 



s (^(1 - X2)mB/V2, br 



1 

"A 



21nJm + ln.^^^^/^) 



ln(6iA) ' "'-ln(62A) 
Sef{t) = s (^XimB/V2,bi^ + s (^X2'mB/V2,b2^ + s (^{1 

+s {x^niB/ v^, 62) + s (^(1 - X3)mB/ V2, 62 



- X2)mB 



/V2,k 



1 

"A 



In 



ln(t/A) 



2 In 



ln(t/A) 



ln(6iA) -ln(62A) 

Sghit) = S {x2mB/ V2, 62) + S [x^VTiBl^Pl, 63 j + s (^(1 - X2)mB/V2, br 

+ s({l- X3)mB/ V2, 63) - 



- ln(t/A) ln(t/A) 
-ln(6iA) -ln(62A)_ 

where the function s{q, b) are defined in the Appendix A of Ref. [2 
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